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Abstract 

We seek analytic formulae for the pricing of extinguishable cap-floor Libor quanto flows which pay 
capped and/or floored foreign currency Libor conditional on survival of a named debt issuer. We assume 
that the foreign currency in which the Libor is paid and the credit default risk (observed for debt 
denominated in the domestic currency) are potentially correlated with each other and with the exchange 
rate between foreign and domestic currency; and further that this exchange rate may jump by a fixed 
relative amount in the event of a credit default event. The stochastic credit intensity model is assumed to 
be Black-Karasinski and the interest rate model to be beta-blend, a family of models which encompasses 
within its scope both the Black-Karasinski and the Hull-White cases. 

The solutions are obtained as perturbation expansions which are valid in the limit of the foreign 
interest rate and the credit default intensity being small. A Green’s function solution is found for the 
governing PDE which is asymptotically valid under this assumption. This is calculated to second order 
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and used to calculate to first order accuracy the price of extinguishable Libor cap-floor quanto flows, in 
particular taking account of the assumed correlations. When the interest rate model is Hull-White, a 
composite expansion is proposed which has the advantage of greater consistency with the known exact 
formula for the conditional bond price in the limit of zero correlation. 


1 Introduction 

We consider an economy with two currencies, foreign and domestic, and a credit default intensity process 
associated either with a corporation based in the foreign economy or a foreign sovereign issuer, the CDS 
market for which operates in the domestic currency. We suppose that the foreign interest rate is stochastic 
but, in the interests of tractability, that the domestic interest rate is deterministic. We also suppose that 
there is a stochastic exchange rate between the two currencies and further that upon default of the named 
issuer there is a decrease in the value of the foreign currency giving rise to a proportional downward jump 
in the exchange rate from foreign to domestic currency. 

A similar model was first considered by Ehlers & Schonbucher (2006), who took the interest rate and 
credit intensity models to be affine, viz. in practice either of Hull & White (1990) or of Cox, Ingersoll & 
Ross (1985) (CIR) type. They did not solve directly the equations they derived but inferred jump levels 
from market data in the light of their modelling approach. Related work was reported by EL-Mohammadi 
(2009) who supposed the credit intensity to be lognormally distributed (with no mean reversion) and used 
his model to price defaultable FX options. A Black & Karasinski (1991) model was first considered in this 
context by Brigo et al. (2015) in their work reviewing the evidence from the comparative USD and EUR 
CDS rates for Italian sovereign debt of an implied FX jump in the wake of the euro crisis. Both the previous 
authors calculated prices using finite difference approaches. In more recent work by Itkin et al. (2017), both 
the interest rates were taken to be stochastic, governed by a CIR model, the foreign rate allowing of a jump 
at default as with the FX rate; the credit model was again taken to be Black-Karasinski. 

A mathematical model effectively identical to that of Brigo et al. (2015) was proposed by Turfus (2016) 
to price contingent convertible (CoCo) bonds, with the equity price taking the role of the FX rate and a 
conversion intensity taking the role of the credit default intensity. Analytic formulae were derived in this 
case based on a perturbation expansion approach. Taking the lead from this work, we have like Itkin et 
al. (2017) sought to extend the modelling approach of Brigo et al. (2015) to include a stochastic (foreign) 
interest rate model, so as to allow cap-floor floating rate coupons to be priced analytically in a defaultable 
quanto framework. However the interest rate model is here assumed to be beta-blend, which approach 
encompasses both the Hull & White (1990) and Black & Karasinski (1991) short rate models. Thus we 
look to price in domestic currency an instrument which pays the foreign Libor rate, possibly capped and/or 
floored, conditional on survival of the named issuer. We will refer to an instrument which pays such cash 
flows as an extinguishable Libor cap-floor quanto (ELCFQ), and in the absence of any cap or floor as an 
ELQ. 

We start in section 2 by defining the underlying processes in terms of suitably chosen auxiliary variables 
and the no-arbitrage conditions they are required to satisfy; we infer therefrom the governing PDE for 
contingent claims. Given the intractability of this PDE, we define in section 3 an asymptotic scaling of 
the underlying variables and associated functions based on the assumed smallness of the two short rates: 
foreign interest rate and credit default intensity. This gives rise to a tractable leading order PDE with the 
“intractable” parts isolated as a perturbation: we can then seek solutions as perturbation series. Rather than 
looking to obtain particular solutions directly, a Green’s function for the full PDE is sought as a perturbation 
expansion, details of which are given in Appendix A. Consistency conditions are derived at this point to 
ensure the model satisfies the required no-arbitrage condition in relation to model-implied prices for risk¬ 
free and risky cash flows. This Green’s function is then used in Appendix B to derive asymptotically valid 
prices for ELCFQ. The main results of the present work are set out in section 4, including some special case 
formulae applicable for the important Hull-White case. Some numerical tests of our formulae against results 
obtained by Monte Carlo simulation are presented in section 5. Conclusions are set out in section 6. 
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2 Stochastic Modelling 

2.1 Definition of Underlying Processes 

Our modelling approach is to represent the foreign interest rate and the credit default intensity as short rate 
diffusion processes and the FX rate as a jump diffusion. Specifically, we suppose the foreign interest rate 
process r t to be governed for t > 0 by a beta-blend short rate model as set out in Turfus (2016), Horvath et 
al. (2017) and Turfus (2017). This encompasses both the Hull & White (1990) and Black & Karasinski (1991) 
short rate models. We suppose the credit default intensity process X t to be governed by a Black-Karasinski 
(lognormal) short rate model. The exchange rate Z t from foreign into domestic currency we take to be given 
by a jump-diffusion process, with a downward jump of a fixed relative amount k < 0 occurring at the time 
of a default of the named issuer. Such jumps are commonly inferred from market data in relation to credit 
issuers which are perceived to carry systemic risk, particularly but not exclusively in the context of emerging 
markets. The diffusive processes are all potentially correlated. We effectively assume no correlation of any 
variables with the domestic interest rate, which justifies our taking it, for computational and notational 
convenience, to be deterministic. 1 Thus we take the domestic interest rate to be given by rd{t), assumed 
known from market data. We have in summary 

Td{t) : domestic interest rate; 

rt : foreign interest rate; 

A t : credit default intensity for debt denominated in domestic currency; 

Z t : exchange rate from foreign to domestic currency. 

We shall find it convenient to represent rt and X t using auxiliary variables Xt and y t , respectively, satisfying 
the following Ornstein-Uhlenbeck processes: 

dx t . = — a r x t dt + a r (t) dWl, (2.1) 

dy t = - uxyt.dt + a\(t ) dW?, (2.2) 

where dW} and dW^ are Brownian motions under the foreign currency and domestic currency equivalent 
martingale measures respectively with 

corr(W t 1 , Wf) = p r \ 

assumed and x 0 = f/ 0 = 0. These auxiliary variables are related to the foreign interest rate r t and the credit 
default intensity Xt, respectively, by 

(1 -P)r t + fdr(t) = (r(t) + (1 - fd)r*(t))£ (^=^) , (2.3) 

X t = (X(t) + X*(t))£(y t ), (2.4) 

with (3 £ [0,1) assumed. Here r(t) is the instantaneous forward rate of foreign interest and A(£) the associated 
credit spread (see (2.15) below), with a r (t) and cr\(t) their respective volatilities and a r and a\ their mean 
reversion rates. Here £(X t ) exp (Xf t — ^[AT] t ) is a stochastic exponential with [X] t the quadratic variation 2 
of a process X t under the requisite measure. The required form of the configurable functions r*(t) and A *{t) 
is determined by calibration of the model to satisfy the no-arbitrage conditions set out below. As can be 
seen, the beta-blend model for the interest rate represents a hybrid between the Hull-White model (/3 f 1) 
and the Black-Karasinski model (f3 = 0). The credit intensity model is Black-Karasinski. 

1 In the absence of correlation between the domestic interest rate and any of the other three stochastic processes, stochasticity 
of the domestic interest rates can be shown to have no impact on the prices we calculate. 

2 Strictly, we will interpret [X]t in this definition as referring to the contribution only from the diffusive part of Xt, ignoring 
that arising from any jump component. 
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Finally we suppose the FX rate to be given by 

dZ 

= (r d (t) -r t - kX t ) dt + a z {t) dWf + kdn t . (2.5) 

where dWf is a Brownian motion and n t a Cox process with intensity A t under the domestic currency 
equivalent martingale measure. We further suppose 

corr(W t 1 , W t 3 ) = p rz , 
corr (W?,W t 3 ) = p Xz . 

It will be convenient to express Z t also through an auxiliary variable, Zt, defined such that 

Z t = F(t)£(z t ), (2.6) 

where 

F{t) := Z 0 e^ (r<i(s) " F(s) " fcA(s)),is , (2.7) 

whence we infer 

dz t = — (r t — r(t) + k(X t — X (t))) dt + <r x {t) dWf + kdn t ■ (2.8) 

with zq = 0. 

We will also need the representation of dx t under the domestic currency equivalent martingale measure. 
This is seen by application of the Girsanov theorem in the standard manner to be given by 

dx t = ~a r x t dt — p rz a r (t)a z {t)dt + oy(t) dWp, (2.9) 

for some new Brownian motion (/IF, 1 related to W.} by 

dWf = dWl — p rz a z (t)dt 


This has solution subject to Xq 


where 


= 0 given by 


x t = -I rz (t)+ [ e s) a r (s)dW^ 
Jo 

I rz (t) ■= Prz [ e~°‘ T ^ t ~ u) a r {u)a z (u)du. 
Jo 


( 2 . 10 ) 

( 2 . 11 ) 


The no-arbitrage conditions 

The formal no-arbitrage constraints which determine the functions r*(t) and A *(t) are as follows. First, by 
considering a risk-free foreign currency cash flow at time t, we deduce 

E f = D(0,t), (2.12) 

for 0 < t < T m , where T m is the longest maturity date for which the model is calibrated, and 

D{t 1 M) = e~ S % ns)ds (2.13) 

is the fi-forward price of the maturity zero coupon bond denominated in foreign currency. Likewise, 
considering a risky domestic currency cash flow, we have 

E d r e -/o(^W+ A ») ds l =B{0,t), (2.14) 

where _ 

B(t i,t 2 ) = e ~ ft? (2.15) 

is the price of a risky domestic currency cash flow. Here E d [-] and [•] represent expectations under the 
domestic and foreign currency martingale measures respectively. We shall assume the bond prices can be 
ascertained at the initial time t = 0 from the market, whence we can view (2.13) and (2.15) as implicitly 
defining f{t) and A (t). 
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2.2 Derivation of Governing PDE 

We consider in the first instance the general problem of pricing a cash security with maturity T whose payoff 
depends on Xt, Vt and zt ■ We introduce the convenient shorthand notation that, for a process X t and 
deterministic function /(.), 

Sx(mX t )) := S(f(t)X t ))\ Xt=x , 

in terms of which we can re-write (2.3) and (2.4) as ry = r(xt,t) and X t = X where 

r(x,t) := j-!—j ^(r(t) + (1 - 0)r*(t))£ & ^ ^ , ( 2 -!6) 

X(y,t) := (X(t) + X*(t))£y(y t ), (2.17) 

Writing the price of the security at time t £ [0,T] as / t T = f(x t ,yt.,Zt,t), we can infer by application of the 
Feynman-Kac theorem to (2.1), (2.2) and (2.8) in the standard manner that the function f(x,y,z,t) satisfies 
the following backward diffusion equation: 

(Jjj. + £ - r d{t) - A(f)^ f(x, y, z , t ) = 0, (2.18) 


where 


d d d 

£■=- a rXgZ - «A y-gZ - (r{x,t) - r(t) + k(X(y,t) - A(f))) — - (A (y,t) - X (£))■ 

+ l ( <r ' (t) ^ + + "-- <f) & + 2P '^ AtMt) idj + J) (l + |) 

+ 2p Xz a x (t)a z {t)-^^j (2.19) 

with in general = P(xt,Zt) for some payoff function P(.). 3 In the absence of closed form solutions to 
(2.18) and guided by the work of Hagan et al. (2015) and Turfus (2017), we propose a perturbation expansion 
approach as follows. 

3 Asymptotic Modelling 

For both short rate models we formally apply a “low rates” assumption. To this end we define small 
parameters 


and 0(1) functions 




r(t) dt , 


£A := 


1 

tt A Pfri 



(3.1) 

(3.2) 


f(t) := e r 1 

erV(f) 

3 It is equally possible to handle a payoff on default within our proposed framework, although we shall not explicitly consider 
that possibility here. 
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and 


\*(t) = e x 1 \*(t). 

Note, we make no formal requirement that either <J r (t) or a\(t) be small (scaled, say, in terms of their 
respective mean reversion rates). But, if they are, this will in general improve the accuracy of the perturbation 
expansions deduced. We further define a new scaled variable Xt and associated scaled volatility a x (t) by 

x t := e~ 0 {x t + I rz {t))e art , 


c(t) = 




(3-5) 

(3-6) 

(3.7) 


Notice that we can also write 

x t := e^e~ art (x t + l xz (0,t)), 
with I xz (.) given by (A.12) below. Likewise we define 

Vt ■= Vte axt , 
v v {t) ■= cr\(t)e axt . 

Here the exponential time scaling is to facilitate removal of the mean reverting drift terms in (2.19). We 
further define new functional forms /(•) and P(-) by: 

f(x,y,Z,t) ■= ft \x t =x,y t =y,zt=zi 

P(x T ,z T ) = P(x T ,z T ), 

where x t and y t are related to xt and y t by (3.5) and (3.7), respectively for 0 < t < r A T. In this notation, 
(2.1) and (2.2) can be re-expressed as 

dx t = a x (t) dWt , (3.8) 

dy t = a y (t) dW 2 (3.9) 

and (2.18) as 


d d 

— + £ - (r d {t) + A (t)) - (e r h(x, t) + e x kg{y , t)) — - e x g{y, t ) ) f(x, y, z,t) = 0 


(3.10) 


where 


: = ^[ a l( t )^+ a l( t )§ y 2+ a z( t )^2+ 2 PrX^x(t)(Xy(t)^ Kr +2prza x (t)(T z (t) 


(p_ 

' dz 2 


d 2 

' dxdy 


d 2 

dxdz 


+ 2p\z<x v (t)a z (t) 


d 2 \ 
dydz J 


is a pure diffusion and 


h(x, t) := h{x, t 1 1 ), 
g(y,t) ■■= g{y,t,t ), 


(3.11) 


(3-12) 

(3.13) 


with 


h(x, t, t\) := 


Fp{t) ■■= 

g(y,t,t i) := 


+ r*{t\))£ x (Fp(ti)(x t — I xz (0, t))) - r(fi), 
\r(t)\P 

(A(ti) + \*(ti))£ y (e axtl yt ) — A(ti), ti>t. 


h>t, 


(3.14) 

(3.15) 

(3.16) 


We seek a Green’s function solution for (3.10) as a joint power series in e r and e x . Details are provided in 
Appendix A. This Green’s function can then be used to generate asymptotic approximations to ELQ and 
ELCFQ products in Appendix B. The main conclusions of this analysis are summarised, in original unsealed 
notation, in section 4 below. 
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4 Main Results 


4.1 Libor Flows 


Theorem 4.1. The price in domestic currency of an Extinguishable Libor Quanto payment at time T of 
tenor-r foreign Libor fixed at time T — t is under the modelling assumptions of section 2 given asymptotically 
in original unsealed variables for 0 < (3 < 1 by 


V, 


ELQ 


F(T)B{0,T) 
D(T — t, T) 


[(1 - D(T - t,T))( 1 - (1 + k)f Xz (T)) - (1 + k)cf rXz {T - r, T)}, 


(4.1) 


with errors = 0(e r (e r + e^)), where 


MT)~J A(v)(e'>‘ M - l)dv, 

MT„T 2 ):= h r * r -» - 1 dudv 

JTi Jo -L — P 


r T i 


rT2 


+ / X(v) / r(u) 

JTi 


_ exp (F 0 (u,v)I r \(v)) - 1 


to 


1 - 


du dv 


Fp{u,v) := 


(1 - ( 0) e -«r(«-») 

Ir(u)!/ 3 


I du, 


I r \{t) ■= Pr\ [ e (“r+a*)(t u la r {u)<T\(u) ( 

Jo 

I\z{t) ■■= p\z [ e~ ax{t ~ u} cr x (u)a z (u) du. 

Jo 

Proof. The derivation of this result is covered in Appendices A and B, in particular at section B.l. 


(4.2) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

□ 


Notice that, in the absence of any correlation involving the credit intensity, the 0.(-) factors disappear 
and the standard result expected under the assumption of deterministic interest rates and credit default 
rates is recovered. The assumption of an exchange rate jump at default continues to give rise to the slight 
modification to the usual definition of the FX forward in (2.7). We infer from this that the 0(e") errors 
propounded above disappear for the special case of t = 0, and the error in (4.1) is in fact 0(e r e x ). 

Although the expression (4.3) is not valid for the Hull-White case ft = 1, it can yet be evaluated in the 
limit as /3 f 1. We have: 


Corollary 4.1. The result in Theorem f.l is also applicable for the case j3 = 1 provided (4.3) is replaced by 
its limiting value 


with 


(T U T 2 ) := l 2 A(n)e-“^ Tl > (e 1 ^ - l) dv f ' 7(m, Ti)I rX (u /\Ti)du 

J T\ J 0 

+ B*(T 2 — Tj) [ *A {v)e~ a ^-^I rX { V ) dv , 

Jo 


(4.7) 


B*{t) 
7 (u,v) 


1 - e - a - T 
a r 

f e -a\(v-u), u < v> 
| e -a r {u-v), u>v5 


(4.8) 

(4.9) 
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Proof. The above results are established by a straightforward application of l’Hopital’s theorem to the 
relevant integrands. □ 


Note that a pure rates-credit correlation effect arises only from the second line of (4.3): a (small) non¬ 
zero impact arises from the first line only if credit intensity and FX are also correlated. For typical market 
parameters and investment grade credits, the PV impact from (4.7) is unlikely to exceed about 1% in relative 
terms, even with the most challenging case of long-dated ly Libor rates. 

Although we do not offer formal justification here, the structure of the Green’s function generator in 
(A.l) suggests that the linear correction involving 4>\ z (-) is in fact the leading order approximations to an 
exponential function. On this basis we propose that a potentially more accurate composite expansion can 
be obtained by completing the exponential as follows: 

Proposition 4.1. A better approxiination than (4.1) should be achieved by using 


Velq 


F(T)B(0,T) 
D(T — r, T) 


e -(i+fc)0A.(r)( 1 _ £,( T _ T)T )) - (1 + k)<fi r \ z (T 



(4.10) 


which has accuracy of the same order. 


4.2 Impact of Caps and Floors 

We consider next the impact of imposing a Libor cap on the value of an ELQ flow. 

Theorem 4.2. The price in domestic currency of a cap at level Kq applied to an Extinguishable Libor 
Quanto payment at time T of tenor-r foreign Libor fixed at time T — t is under the modelling assumptions 
of section 2 given asymptotically in original unsealed variables for 0 < /3 < 1 by 


Vcaplet 


F{T)B{0,T) 
D(T — t, T) 


(1 - (1 + k)f Xz (T)) (1 - (1 + K C S(T - t, T))D(T - r, T^N^d^x*,T - r)) 


<' T N(-d 2 (x*,u,T-T))-N(-di(x*,T-T)) 

r[u) ----- au 


IT-t 


i -p 


- (1 + k)ip,(x*,T-T,T) 


(4.11) 


with errors = 0(e r (e r + e\)), where S(-) represents the day count fraction calculated according to the relevant 
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convention, 


^ Xz (x,T u T 2 ) := p \{v) - l) f(u) 

J T\ J 0 


exp(.F/3(w, u A Xi)e uATl) I r \(u A Ti))A/’(±rf 3 (f*, u, i>, Ti)) - AT(±d 2 (x*, u, Ti)) 

1 -I 

_ exp(F / 3 (u, v)I r x(v))N(±d 3 (x*,u, v, Ti)) - N(±d 2 {x*, u, Ti)) 


dudv 


rTi 


pT2 


+ / A(v) / r(u) 
Jo J Ti 


1 -, 


I r {t) := f e 2oir ^ U ^a 2 (u)du, 
Jo 


di(x,Ti) := 


d 2 (x,u,Ti) := di(x) - 


d 3 (x,u,v,T 1 ) := d 2 (x,u,Ti) - 


s/Wi)' 

Fp(u + Ti, 2(u A Ti))I r {u A Ti) 


V^r(?T) 

Fp(u,u Av A Ti)e~ a ^^ v ~ uAvATl ^ I r \(u A v A T x ) 


(in du, 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


/T-t 


and x* satisfies 

r T 

f(s) exp(Fp(s,T—T)x* — ?iFp(s)I r (T—T)) ds = (1—P)(1 +KcS(T—t,T))D(T—t,T)+P—D(T—t,T)~ 1 . 

(4.17) 

Proof. The derivation of this result is covered in Appendices A and B, in particular at section B.2. □ 

Note that the <j>\ z {-) term can again if wished be replaced with an exponential. The Hull-White result is 
again obtained by letting ft f 1: 

Corollary 4.2. The result corresponding to that in Theorem f.2 for the case /3 = 1 is 


V r 


hw F(T)B(0, T) 


caplet 


D(T — t, T) L 


(1 ~ (1 + k)fi Xz (T)) (J1 - (1 + K C S(T - r, T))D(T - r, T))N(—d\{x*,T - r)) 

+ B* (r) \JI r {T - T)N'(—di (x*,T — r))) - (1 + k)^ w (x*,T- t,T) 

(4.18) 


where 

ip^ W (x,T 1: T 2 ) := <j>™{T u T 2 ) (N{±d 1 {x,T 1 ))=f B* (r)^ I r (T - T )N , {±d 1 (x,T 1 ))) , (4.19) 

and the required value of x* is 

,* (1 + K C S(T — t,T))D(T — t,T) — 1 

* B*(t) 

Proof. This is a straightforward application of l’Hopital’s theorem. □ 

Alternatively, following the logic of Proposition 4.1 above and with an eye on the known exact Hull-White 
solution, we suggest: 
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Proposition 4.2. A better approximation than (4.18) for the case 0 = 1 will be achieved by using the 
composite expansion 


yHW 
* caplet 


F{T)B( 0,T) 
D(T-t, T) 


e -(i+k)4>xz(T) 


(N(-dc) - (1 + K c 6(T 


T,T))D(T-r,T)N(-d +)) 


(l + k)<j>™(T-T,T)N(-dc) 


(4.20) 


with 

ln(£>(r - r,r)( 1 + KqS(T - t,T))) ± \B*{rfl r (T - r) 
C ' B*(r)\J I r (T — t) 


(4.21) 


The expression (4.20) can be shown to agree with (4.18) to 0(e r e\), but has the desirable property of 
matching the known exact solution in the limit of zero correlation. Horvath et, al. (2017) offer evidence that 
the difference between results arising from the two approaches is typically very small in practice. In fact, 
as with the result in Theorem 4.1, because the known exact result is recovered when the </>.(•) factors go to 
zero, we infer that errors in (4.20) are likewise 0(e r e^). 

It is also worth noting that, comparing (4.20) with the result arising in the absence of an assumed jump 
at default and stochastic credit, the impact of the jump at leading order is to multiply the effective credit 
intensity A(f) by a factor of 1 + fc, while the impact of FX-credit correlation is to multiply by a further factor 
of e Ixz< ' t ' ) . This should be equally true for all values of 0. The impact of rates-credit correlation is of course 
not so readily captured. 

By a similar means or, more straightforwardly, by applying the put-call parity principle to (4.11), we see 
that the corresponding floorlet price for 0 < 1 with floor at Kp is, in original unsealed notation, given as 
follows: 


Corollary 4.3. The floor price corresponding to the situation in Theorem f.2 is 


Vfloorlet 


F(T)B(0, T) 
D(T — r, T) 


(1 - (1 + k)4>\ z (T)) ^((1 + K F 5(T - r, T))D(T — r, T) — 1 )N(d 1 (x*,T - r)) 

[ T ^N(d 2 {x*,u,T — r)) — N(di(x*,T — r)) A 
— / r[u) ----- au 


/ T-7 


1-/3 


+ (1 + *Wa x {x*,T-t,T) 


(4.22) 


with errors = 0(e r (e r + e^))- 

Note that the (j>\ z (-) term can again if wished be replaced with an exponential. Similarly we have: 
Proposition 4.3. The floor price approximation corresponding to (4.20) for the case 0 = 1 is 


V t 


HW _ F(T)B{0,T) 
D(T — t, T) L' 


floorlet 


e -(!+fc)<K,(T) ^ + Kf§ ( T _ T) T))D(T - T, T)N(dp) - N(dp)j 

+ (1 + ( T - r, T)N(dp)} (4.23) 


with 


d$ := 


ln(D(T - t,T)(1 + K f 6(T - t,T))) = ±\B*(r) 2 I r (T - t) 

B *(. T )\/Ir{T - t) 


(4.24) 


Finally we observe that the value of an ELCFQ capped at Kc and floored at Kp < Kq can be priced as 
a combination of an ELQ with a caplet and a floorlet, viz. 


VeLCFQ = VeLQ — Vcapiet + ^floorlet- 


(4.25) 
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5 Numerical Tests 


We test the formula (4.25) with (4.20) and (4.23) for an ECLFQ paying 3m Libor over five years. The foreign 
yield curve is assumed flat at 2.5%. A floor is imposed on Libor payments at 1.5% and various cap levels are 
considered. We take the Hull-White (normal) volatility to be ay = 1% with a mean reversion rate ay = 0.25, 
the 5y credit spread to be 40bp with a (lognormal) volatility of a\ = 0.65 and a mean reversion ot\ = 0.3, 
and the FX volatility to be oy = 0.25. The correlation between FX and credit is taken to be p\ z = 0.4 
with a supposed downward jump of k = 0.091 in the FX rate at default. Comparisons between our analytic 
formulae and numerically computed prices based on Monte Carlo simulation are shown in Fig. 1. As can 
be seen the agreement is excellent and indeed comparable to the numerical error achieved with 500,000 MC 
simulations. 
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Figure 1: ELCFQ prices for various cap levels using Hull-White interest rate model. 
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Figure 2: ELCFQ prices for various p\ z levels using Hull-White interest rate model. 

A test is also done showing the impact of imposing a correlation p\ z between the credit intensity and the 
FX rate. The details of the calculation are as above except that no cap is imposed, only a floor. Results are 
illustrated in Fig. 2. The impact of the correlation appears to be well captured by our analytical model, to 
within the numerical error associated with the Monte Carlo simulation. 


6 Conclusions 

We have computed approximate formulae for the pricing in domestic currency of extinguishable cap-floor 
Libor quanto flows which pay capped and/or floored foreign currency Libor conditional on survival of a 
named debt issuer. We assumed that the foreign currency in which the Libor is paid and the credit default 
risk (observed for debt denominated in the domestic currency) are potentially correlated with the exchange 
rate between foreign and domestic currency. The domestic interest rate was taken for convenience to be 
deterministic. Solutions were obtained as perturbation expansions valid in the limit of the foreign interest 
rate and the credit default intensity being small. A Green’s function solution to the governing PDE was 
derived accurate to second order in both parameters although in practice we only needed to use the first 
order terms of the Green’s function to obtain what we believe will be adequate accuracy for many if not 
most circumstances. 

In the absence of any cap or floor, the value of the Libor flow is given by (4.1) or (4.10). In the event 
that the flow is capped at a level Kc and/or floored at a level Kp , the value must be adjusted in accordance 
with (4.25). The caplet and floorlet formulae needed in this context are given by (4.11) and (4.22) above, 
respectively, for /? £ [0,1). Variants of the formulae suitable for use in a Hull-White context (/3 = 1), 
which can be considered a limiting case of the aforementioned formulae, are given by (4.20) and (4.23). 
Comparisons of these formulae with results obtained by Monte Carlo simulation show them to be highly 
accurate. 
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A Green’s Function 


A.l Leading Order Solution 

We follow Hagan et al. (2015) in observing that the operator 


f v ( d 

U(t;v)=ex p J \£(u) - e r h(x,u) — - e x g(y,u) 


1 + k^~ 
oz 


du 


(A.l) 


is a formal solution of 


( e rh{x , t) + e A kg(y, t)) - e\g(y, t)^ f(x, y, z,t) = 0 (A.2) 

for t < v subject to U(v\v) = /. The Green’s function for this problem is then the integral kernel of U(t\v ), 
viz. 

G{x,y,z,t;£,r},C,v) = U(t;v)(x,y,z;^,‘n,Q. (A.3) 

The Green’s function for (2.18) will then be given straightforwardly by 

G*(x,y,z,t;£,r},v) = B(t,v)G(x,y,z,t;^,r], (,v). (A.4) 


We concentrate henceforth on the calculation of G(-). Again following Hagan et al. (2015) we consider first 
the limiting problem with e r = e\ = 0 as the leading order term for a perturbation expansion. The required 
solution can be written 

/ V 

C{u)du. (A.5) 

We will in all cases be interested in so-called free-boundary Green’s function solutions which tend to zero as 
x,y —> ±oo. The Green’s function solution subject to these conditions is well known. It is given by: 

g3 

Go,o(x,y,z,t-,€,r),(;,v) = g^Q^N 3 ((;-x,r)-y,C- z;R{t,v)), t<v (A.6) 


where Ns(x, y, z; R{t, v)) is a trivariate Gaussian probability distribution function with mean 0 and covariance 
matrix 


I xy (t,v ) I xz {t,V) \ 

I X y(t,v) Iy(t,v) Iyz(t,v) 
Ixz{t,V) Iyz(t,v) Iz(t,V) J 


(A.7) 


with 


t-2) 

I X y (fli t‘2 ) 
Ixz(t l,t-2) 
dyz (^1i^ 2 ) 



al(u)du, 
a 2 y (u) du, 
<J 2 (u) du, 


ft 2 

<Jx{u)tJy{u) du, 

Jtx 


pt2 


/ 

ffi(ti)ff 2 («) du, 

Jt! 


/**2 


/ 

a y (u)a z (u) du. 

Jt, 
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(A.8) 

(A.9) 

(A.10) 

(ATI) 

(A.12) 

(A.13) 



Before deriving higher order terms in our proposed Green’s function expansion, we introduce some additional 
notation. First, following Turfus (2017), we propose that the form of the configurable function f*(t) required 
to satisfy (2.12) is 

OO 

r*(f) = £4r*(f). (A.14) 

i = 1 

Considering next the form of the expansion for A *(t) required to satisfy (2.14), we note that Turfus (2017) 
considered the case where rates were stochastic. Here domestic rates are assumed deterministic so the 
calculation is on this occasion simpler (independent of the interest rate expansion). We are led to expect 

OO 

V(t) = '£&](*)■ (A-!5) 

i=i 


On this basis we can now expand h(x,t,ti) as 


h(x, t, ti) = ^2 *i) - r(ti) (A.16) 

i =0 

where 


and similarly 


where 


ho(x,t,ti) := r(ti)£ x (Fp(ti)(x t - 4*(M))) 
hi(x,t,ti) := r*(ti)£ x (Fp(ti)(x t - l xz (0,t))), i > 0. 


9 {y,t,t i) ='^2e J x g j (x,t,t 1 ) - A(G) 
i—o 


9o(y,t,ti) := A(h)£y(e axtl y t ) 

9 j(y,t,t i) := X*(t 1 )£ y (e~ axtl y t ), j > 0. 


(A.17) 


Note that the required quadratic variations are given by \x\t = 4(0, t) and [y\ t = 4(0, t). We further define 


hi(x,t) := hi(x,t,t), i> 0, 

9j(y,t) ■= 9 j(y,t,t), j> o. 


(A.18) 

(A.19) 


A.2 Asymptotic Expansion 

Taking into consideration the structure of (3.10) and following Kato (1995, see in particular Chapter IX, 
§2), we propose that a formal asymptotic expansion of U(t\ v) is possible as 

OO OO 

U(t;v) = ^2^2e l r e{Uij{t;v) (A.20) 

i—0 j—0 

where each of the is 0(1). By a process of induction we infer that the following recurrence relation 

must hold: 


77 / ,'j (t. v ) 


W 0 ,o (t;u) 


h(x,u)—U t - ld (u; v)l i>0 


g(y,u) 


1 + k^~ 

az 


Uij-i(u;v)l J>0 du 


i,j > 0. 

(A.21) 
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Substituting this expression into (A.20) and and grouping like powers of e r and c\. we can write 


U(t-v) = ^2^2e l r e{U* J (t;v), 

2 — 0 j= 0 

where in particular we have 

Uo,o( t ’ v ) = w 0,0 

[ v a 

U*,a(t\ v ) = - J U 0 ,o{t;u)(ho(x,u) - f(u))—U 0 fi(u;v)du, 

Ko,i(t',v) = - W 0i0 (t; u)(go{y, u) - A (u)) ^1 + k-^j U 0 , 0 {u; v)du, 

[ v ( d d \ 

G 2fi {t;v) = - J U 0 ,o(t;u) I (h 0 (x, u) - r(u))— U lfi (u; v) + hi(x, u) — U Q , 0 (u; v) \ du, 


d 


U*,i(t\ v ) = ~ I Uo,o(t;u) ( (h 0 {x,u) - r(u))—U 0 ^(u\v) + ( g 0 (y,u) - A(it)) ( 1 + fc— ) U lfi (u;v) ) du 


d 


d 


U^ 2 (t\v) = - / Wo,o(t;u) ( {go(y,u) - X(u)) ( 1 + fc— ) Uo,i(u;v) + gi(y,u) ( 1 + k— ) U 0 , 0 (u;v) ) du. 


Taking Gjj(.) to be the integral kernel of Id* •, we deduce: 


a 


a^ 


oo oo 


G*(x, y, z, t \£, rj, (, v) = B(t, v) ^ ^ e % r ^ x G itj {x, y, z, t; £, rj, C, v). 

2 — 0 J=0 


(A.22) 


We seek to obtain explicitly the form of the required integral kernels by a recursive process. We note in this 
regard that an analogous calculation was performed by Pagliarani et al. (2011) who addressed the problem 
of extending the Black-Scholes Green’s function to address local volatility. At first order we obtain 

pV pOO poo p oo 

Gi,o{x,y,z,t;^,g,C,v) = - / / / / G 0 ,o(x,y,z,t;x 1 ,y 1 ,z 1 ,t 1 )(h 0 (x 1 ,t 1 )-r(t-i)) 

J t J—oo J— oo J— oo 


dz\ 


Gofl{x\,yi,z\,t\-, £, rj, (, v ) dx^dy^dz^dti 


r a 

= -j (ho{x,t,t 1 )Ml tl -r(t 1 ))—Go,o(x,y,z,t;€,r),(;,v)dt 1 (A.23) 


and 

Go,i{x,y,z,t;^,r],C,v) = - 


poo poo poo 


Go,o(x,y,z,t;xi,yi,zi,ti)(go(yi,ti) - A(ti)) 


i «/ — oo J — oo J — oo 


a A 

1 + k—J Go,o(xi,yi,zi,ti \£, T), C, v) dxidyidzidti 
= - [ {go(y,t,t 1 )Ml t -X(t 1 ))(l + k^\G 0 ,o{x,y,z,t;£,r),C,v)dt 1 , (A.24) 


a* 


where we have defined the shift operators 


■Mt ut 2 G 0 ,o(x, y, z, t ; £, g, C, v) := G 0 , 0 (x + Fp{t 2 )I x (t 1 ,t 2 ), y, z, t; £, g, <, i>) (A.25) 

y, z, t; £, g, (, v) := G 0 , 0 (x, y + e~ axt 2 I y (t 1 ,t 2 ), z + I yz {t ly t 2 ), t; £, 77 , C, v) (A.26) 
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Similarly at second order we obtain after some manipulations 

C v pt2 


n t-2 

- r(t 1 ))(h 0 {x,t,t 2 )M x tlM - r{t 2 )) 


-qZ2 G o,o(x, y , 2 :, t ; £, 77 , C, v) dt 1 dt 2 


+ J (h 0 (x,t,t 2 ) J {h 0 {x,t,ti) (exp(F / 3 (ti)F i a(t 2 ) 4 ( 0 ,ti)) - 1 ) 


a 2 


—G 0 ,o{x, y, 2 , t ; £, 77 , C, v) dGdt 2 

r v d 

-J h 1 {x,t,t 1 )Ml tl -^G 0 fi(x,y,z,t-,t,r],(,v)dt 1 , (A.27) 


(yo(y,Mi)A4 t , tl - Kti)){go(y,t,t 2 )M tl ,t 2 - A(f 2 )) 

a A 2 


1 + fc— ) G 0 ,o(z, y, t\ £, 77 , C, «) dtidt -2 


+ J (go{y,t,t 2 ) J (g 0 (y,t,ti) (exp (e ax{tl+t2) Iy(0,t 1 ) S j - 1 ) 

( 1 + fc — j ■Mt,t 1 M tu t 2 Go t o(x, y, z, t ; £, 77 , C, «) dGdt 2 
- £ 9 i(y,t,ti)Mt,t 1 (l + k-^jGo t o(x,y,z,t;€,T),(;,v)dt 1 (A.28) 


and 


Gi,i{x,y,t-,t;,r),v) = J J {ho{x,t,t{)M^ tx - r(U))(g 0 {y,t,t 2 )M v tlM - A(t 2 )) 

d \ d 

1 + ^ 7 ^ ) —G 0 ,o(x : y,z,t-^,ij,C,x)dtidt 2 


dz J dz 


v r t 2 


(. go(y,t,ti)Ml tl — A(ti))(/z 0 (a;, t 2 )A^f 1>42 - r(f 2 )) 


9 A 9 

1 + / ]g z Go ^ x ' y ' z ' t '^^ 1 '^ v ^ dtidt ‘ 2 - 


+ [ h 0 (x,t,t 2 ) [ go{y,t,ti) (exp (e aAtl -F>(t 2 )/ p ( 0 , G)) - l) 

/ 9 \ 9 

( 1 + *^7 J ^Ml tl M% u t 2 Go,o{x, y, z , f; £, 77 , C, v) dhdt 2 

+ f go(x,t,t 2 ) f h 0 (x,t,ti) (exp(e~ axt2 F 0 (t 1 )I p (O,t 1 )) - 1) 

Jt Jt 

1 + k §~ S ) ^ M hi MV tite G 0fi(x,y,z,t\£,'qX,v)dtidt 2 . (A.29) 


Substituting (A.23)-(A.29) into (A. 22 ) gives the required Green’s function to second order accuracy. 

As noted by Pagliarani et al. (2011), since the differential operators appearing in the Green’s function 
apply only to 2 and not to £, 77 or £, they do not influence the applicability of the Green’s function to the 
requisite payoff function, which consequently has no smoothness conditions imposed other than integrability, 


16 



say by its being piecewise continuous and exponentially bounded. As in the case of Pagliarani et al. (2011), 
it is if wished possible to apply the zero-order Green’s function to the payoff and then apply the differential 
operator in 2 (and the integrations over time) thereafter. However, since the payoff functions we shall be 
considering here themselves rely on application of the Green’s function to a pre-specified order of accuracy, 
this option is less relevant for us. 


A.3 Calibration 


It remains to calibrate our model consistent with the 110 -arbitrage conditions (2.12) and (2.14). This is 
achieved by considering in the former case the consistent pricing in domestic currency of a risk-free foreign 
currency cash flow, and in the latter case by consideration of a risky cash flow, as we shall now show. 

We seek to calibrate the foreign interest rate model to ensure that (2.12) is satisfied. The necessary 
calculation for a risk-free cash flow in our model is equivalent to that performed by Ttirfus (2017) except 
here we express the results under the transformed variables for the domestic currency measure. We find the 
PV (in foreign currency) at time t is given by fj = X T (x t — I xz (Q,t),t) where 

X T (x, t) ~ D{t, T) (l + e r F lfi (x, t) + e 2 r F 2fi {x, t)) (A.30) 

with 0(e x ) errors, and 





r(G)(£: 2 .(F ) 3(ti)x t ) - 1) dti 


F 2fi {x,t) := \F 2 0 {x,t) + J £ x {Fp(t 2 )x t ) 


f(t 2 ) / f(ti)£ x {Fp[ti)x t ) (exp {Fp{ti)Fp{t 2 )I x (QM)) ~ 1 )dt 1 -r*(f 2 ) dt 2 . 


Of interest to us here is the conclusion that, for the expression in (A.30) to satisfy (2.12) above, we require 
that .£ 1 , 0 ( 0 ,0) = F 2 , 0 (0,0) = 0, so should choose 


r*(f)=r(f) [ r(u)(exp(Fp(u)Fp(t)I x (0,u))-l)du. (A.31) 

Jo 

Higher terms will not be needed for the purposes of the present calculations wherein we truncate the Green’s 
function at second order. 


Pricing of risky domestic cash flow 

The calibration of the credit intensity model proceeds by considering a risky domestic currency cash flow 
and, because domestic interest rates are assumed deterministic, is analogous to that set out above. For 
details, the reader can consult Turfus (2017). We state here only the conclusion, namely that satisfaction of 
(2.14) requires us to choose 

A*(f) = A (t) J A (u) (exp (e- aA ( t+ “)/ y (0, — l^j du, (A.32) 

Note the form of this is exactly analogous to (A.31) with /3 = 0. This completes the calibration of our model 
to second order. 


B Asymptotic Pricing 

B.l Extinguishable Libor Flow 

We consider first the pricing of an ELCFQ without the cap/floor feature, denoting this Velq(*(, Ut, Zt, t). 
This can be modified to obtain the requisite ELCFQ price by subtracting the value of a caplet and adding 
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that of a floor, as we shall see below. We shall throughout make use of a first order Green’s function 
approximation, i.e. using 

G(-) ~ G 0 ,o(-) + trGifi(-) + e\Go t \(-). 

To the same order of approximation we can ignore the second order term 4.0 (•) in (A.30) above. 

The payoff at time T for a payment of tenor- r foreign currency Libor fixed at time T — r can be written 
in domestic currency as 

Pelq(x T -t =x,zt = z) = F(T)£ z (z t )(X t (x - 4,(0, T-t),T- t)- 1 - 1), (B.l) 

with X T (■) defined by (A.30) above. 4 

Viewed from the perspective of time T — r, the payoff (B.l) is independent of both x and y, so its PV 
conditional on (xt-t,Ut-t, zt~ t ) = ( x , y , z) is seen by straightforward application of our first order Green’s 
function to be 


Velq(x, y,z,T - t) = F{T)£ z {z T -t) 


B(T — r, T) 
D(T — t, T) 


^l-e A (l + fc) J (g 0 (y,T -r,v)ex p (e axv I yz (T - t,v)) - A(v)) dvj 
^1 - D{T - r, T) + e r J ( h 0 (x, T — t,u) - f(u)) du'j 


+ 0(e r (e r + e A )), 


where the 0(4) errors arise on account of us having chosen to ignore second order terms in the payoff from 
the outset. We can now treat this expression as the payoff at time T — r for 1 elq(£> y, z, t) for t < T — r. 
We are interested in particular in the price at t = 0. Applying our first order Green’s function to this final 
value problem we obtain 


r r r 

Ve LQ (* = 0)~ 44^(0, T-t) / // G ili (0,0,0,0;^» ? ,C ) r-r)Vfe E Q(^»7,C ) r-r)ded»7dC 

i,j>0 ^ ^ 

~ F oIt B -tT) ( (1 " D{T - T ’ T))(1 “ £a(1 + k)4>vz(T)) - e r e A (l + k)^ xyz (T - r, T)) (B.2) 
with errors = 0(e r (e r + e A )), where 


4>yz(T) 


fixyziTl, 4) 


A(u)) (exp (e axv I yz {0,v)) - l) dv, 
X(v) (exp (e~ axV I yz ( 0,4) - bf 


+ 


t 2 


X(v) (exp (e axV I yz {0,v)) - l) 


(B.3) 

f(u) (exp (Fp{u)e~ axV I xy {0,Ti)) - l) dudv 
■Ti 

r(u) (exp (Fp(u)e- axV I xy (0, u)) — l) dudv 


Ti 


Jo 



rT2 

X( v ) / f(u) (exp (Fp(u) 
Jn 


e axv I xy (0, v)) — l) dudv. 


(B.4) 


4 We ignore here any possible basis between the Libor curve and that for the risk-free rate captured by rt . However, it is 
not difficult to deal with such a basis by using a discount curve which represents the foreign currency tenor-r Libor rate rather 
than the risk-free rate, since for discounting purposes we always use the domestic currency curve (and we suppose there are no 
domestic currency Libor payments). 


18 



Here the first line in the expression for <j) xyz (-) arises from Go,o> the second from Gpo and the third from 
Go,i- This can be rewritten more compactly as 

</>xyz{Ti,T 2 ) := [ X(v) (exp (e~ axV l yz (0,v)) - l) [ r(u) (exp (f' l s(tt)e““ x ’'4 ! ,(0,«A ?i)) — l) dudv 
Jt i Jo 

+ ( \{v) f r(u) (exp (Fp(u)e- axv I xy (0, i>)) — l) dudv. 

JO Jti 

(B.5) 

Expressions in original unsealed notation for the PV of Extinguishable Libor Flows are presented in section 
4.1 above. 

B.2 Extinguishable Caplet 

To calculate the impact of the Libor calculated above being capped at Kq, we must compute the value of a 
caplet which pays the amount by which the Libor exceeds Kq, if this is positive, conditional on survival of 
the named issuer, multiplied by the relevant day count fraction. Thus we have a payoff at time T given by 

-Pcaplet^T—T =X,ZT = z) = F{T)£ z (z T )(X T (x - I xz ( 0, T — t),T — t)' 1 - K~ 1 )l X T^-i xx ( 0t T-r) > T-r)<K, 

(B.6) 

where 

1 

K := 1 + K c 5{T - r, T) ' 

Proceeding as previously with first order representations, we deduce our first order solution as of time T — r 
to be 


E cap let(z, V, Z,T — t) = F(T)£~(z T - t ) ^ 1x>x*+I^(0,T-t) 


^1 - e A (l + k) J (g 0 {y, T-t,v) exp (e axv I yz (T - r, v)) - A(t 
^1 — k~ 1 D(T — t,T) + e r J ( ho(x, T — t,u) — f(w)) dv^j 


+ 0(e r (e r + e A )), 
where x* is the unique solution of 


[ (h 0 {x + I xz (0,T-T),T-T,t 1 )-r(t 1 ))dt 1 = K 1 D(T - t,T) - 1. 

Jt-t 
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Applying our first order Green’s function with the T — r solution as payoff, we obtain the following expression 
for the PV as of t = 0: 


i+j< 1 


V cap iet(i = 0)~ V cU{B(0,T-t) Gi,, (0,0,0,0;C, 77 , C,r-r) Caplet (Z,ri,C,T-T)dSdT,dC 
i,j> 0 J J JR 3 


F(T)B(0, T) 
D(T — r, T) 


(1 - c A (l + k)(j> yz (T)) ^(1 - k~ l D{T - t, T))N(—di(x*,T — r)) 

+ e r J ?(u) ^N(—d 2 (x*,u, T - r)) — N(—d 3 (x*,T — 


- e r e A (l + k)ip xyz {x*,T - r, T) 


(B.7) 


with 0(e r (e r + e|)) error, where 
rT2 


^xyz(x*,T 1 ,T 2 ) = [ X(v) (exp (e u)) - l) [ r(u) 

J T\ Jo 


(exp (Fp(u)e axv I xy ( 0, u A Ti)) N(—d 3 (x*,u, v , Ti)) — N(—d 2 (x*,u, Tf))) drtdu 




pT2 


+ / A(u) / f(u) 
J 0 J Ti 


(exp ( Fp(u)e axV I xy (0,v)) N(-d 3 (x*,u,v,T 1 )) - N(-d 2 (x*,u,T 1 ))) dudv 


(B.8) 


and 


44,4) := 


d 2 (x,u,T 1 ) := 44) - 


v%(o ,T,y 

F/3 {u)I x (0, u A Ti) 


d 3 (x,u,v,Ti) := d 2 (x,u,Tx) - 


44(0 ,Ti) 

Fp(u)e~ axv I X y(0, «At)ATi) 


44(0,4) 


(B.9) 

(B.10) 

(B.ll) 


Note that by construction 

lim ipxyz (x*,T 1 ,T 2 ) = <t> xyz (4, T 2 ) (B.12) 

a; —>• —00 

Analogous results can be derived for floorlets by appeal to the put-call parity theorem. Expressions in 
original unsealed notation for both caplets and floorlets are presented in section 4.2 above. 
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